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Abstract 

Variational Monte Carlo calculations of the quasielectron and quasihole ex- 
citation energies in the fractional quantum Hall effect have been carried out 
at filling fractions v = 1/3, 1/5, and 1/7. For the quasielectron both the trial 
wave function originally proposed by Laughlin and the composite fermion 
wave function proposed by Jain have been used. We find that for long-range 
Coulomb interactions the results obtained using these two wave functions are 
essentially the same, though the energy gap obtained using the composite 
fermion quasielectron is slightly smaller, and closer to extrapolated exact- 
diagonalization results. 
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I. INTRODUCTION 



Shortly after the discovery of the fractional quantum Hall effect (fractional QHE) 1 1 ] 
Laughlin M introduced a set of trial wave functions describing the 'parent' quantum Hall 
states occurring at Landau level filling fraction v = 1/q where q is an odd integer. In addition 
to the ground state, Laughlin introduced trial wave functions describing fractionally charged 
quasielectron (e/q) and quasihole (—e/q) excitations. From the very beginning it was clear 
that the wave function for the quasihole, with its simple Jastrow form, was more natural 
than the wave function for the quasielectron, which contains explicit derivatives with respect 
to electron coordinates. This difference is reflected, for example, in the fact that while there 
exists a Hamiltonian for which the Laughlin ground state and quasihole wave functions are 
exact (zero energy) eigenstates || there exists no such simple Hamiltonian for which the 
quasielectron wave function is also an exact eigenstate. 

According to the composite fermion theory, proposed by Jain |4j], the fractional QHE 
corresponds to an integer QHE of composite fermions - electrons bound to an even number 
of statistical flux quanta. This identification leads to a procedure for constructing fractional 
QHE trial states by first constructing integer QHE states, then multiplying by a Jastrow 
factor which binds an even number of vortices to each electron, and finally projecting the 
resulting state onto the lowest Landau level. The wave function obtained using this proce- 
dure for the parent quantum Hall state is identical to Laughlin's ground state, and the same 
is true for the quasihole wave function. However, the composite fermion quasielectron wave 
function is not identical to the one proposed by Laughlin. 

To date, the best estimate of the energy gap for creating a quasielectron-quasihole pair 
with infinite separation at v — 1/3 computed using Laughlin's trial states was obtained by 
Morf and Halperin, using the disk || and spherical |J geometries, with the result ~ 
0.092 ± 0.004e 2 /e/o where 1$ = yJhc/eB is the magnetic length. This may be compared 
with the result of Bonesteel using the composite fermion quasielectron wave function of 
A C f — 0.106 ± 0.003e 2 /eZ - It therefore appears that the energy gap computed using the 
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Laughlin quasielectron is over %10 lower than that obtained using the composite fermion 
theory, and so one might conclude that the Laughlin quasielectron wave function provides a 
(slightly) better description of the physical quasielectron state. More recently, Girlich and 
Hellmund || have shown that for the truncated (Vi) pseudopotential interaction introduced 
by Haldane ||, the interaction for which the Laughlin ground state and quasihole state are 
exact eigenstates, the Laughlin quasielectron has an energy which is 18 % higher than the 
composite fermion quasielectron. These authors go on to speculate that the same would be 
true for the Coulomb interaction, and this has motivated us to reexamine the calculations 
of Morf and Halperin. 

Taking advantage of the availability of significantly faster computers, and performing 
a better extrapolation of finite size results to the thermodynamic limit, we find that the 
v = 1/3 energy gap computed using the Laughlin quasielectron extrapolates to A L ~ 
0.110(2)e 2 /e/o, a significantly higher result than previously reported. This result is higher 
than Acf and so is consistent with Girlich and Hellmund ||, though we find that for the 
Coulomb interaction the difference between the two energy gaps is quite small (less than 5 %) 
indicating that both wave functions provide adequate descriptions of the true quasielectron. 
We have also performed what we believe are the first calculations of the energy gap using 
Laughlin's quasielectron wave functions for filling fractions v — 1/5 and 1/7. Comparing 
these energies to the corresponding composite fermion energies we find the same result — the 
composite fermion energy gap is consistently smaller than the corresponding Laughlin energy 
gap. Comparing these results to the extrapolated exact diagonalization results of Fano et al. 
we find that the composite fermion energy gaps are also consistently closer to the 'exact' 
results. However, as for v = 1/3, the differences are slight, and the main conclusion is that 
both wave functions provide an adequate description of the physical quasielectron. 

This paper is organized as follows. In Sec. II we review the formulation of the two 
dimensional electron gas on the Haldane sphere and introduce a procedure for projecting 
wave functions on the sphere into the lowest Landau level. In Sec. Ill we review both the 
composite fermion construction of the quasielectron state, which requires the projection 
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developed in Sec. II, as well as the Laughlin quasielectron wave function. Finally, Sec. IV 
contains a summary of our results. 



II. THE TWO DIMENSIONAL ELECTRON GAS ON THE HALDANE SPHERE 

A. One Particle and Landau Level Projection 

We begin by reviewing the problem of one spin polarized electron confined to move on the 
surface of a sphere of radius R with a magnetic monopole at its center. Following convention 
we let 

denote the number of flux quanta piercing the surface of the sphere due to the monopole. 
The Hamiltonian describing this particle is then 



IA 
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T = ^^- (1) 
where A = r x (— iKV + eA(r)) and V x A = Br on the surface of the sphere. In what 



follows we work in the Wu-Yang gauge |10| for which 



hS 1 — cos 9 , . 

A = (2) 

eR sm9 V ; 

and use the complex stereographic coordinate z = x + iy = cot |e 2 ^. The eigenf unctions of 



T in this gauge are the monopole harmonics |10 



Y Slm = M Slm f j z S + m p S+m,S-m ( _J^_ ) , 3 , 



Or 


z 






z 





where 



m ~( 2l+1 v-^+sy- V /2 u) 

Msim - \ 4vr (/-m)!(/ + m)! ) ' (4) 

P° ,/3 is a Jacobi polynomial, / = S, S+l, S+2, ■ ■ -, and for a given l,m — —I, — l+l, • • • , l—l, I. 
If we let n — S — I then the energies of these states are 

4 



E n — hit) ( n 



1 n(n + 1; 



2S 



n = 0,1,2, 



(5) 



and n is the spherical Landau level index. 

It will be necessary in what follows to project general wave functions onto the lowest 
Landau level (n = 0) Hilbert space. Following Girvin and Jach [ I I | we now introduce a 
general procedure for performing such a projection on the sphere. First note that the lowest 
Landau level wave functions in the Wu-Yang gauge are 

\\ 1/2 , , >. 



/ / 

2S + 1 



V 



47T 



2S 
S + M 



1+ z 



2 *> 



p= 0, 2S + 1. 



The Hilbert space of lowest Landau level wave functions on the sphere then corresponds to 
wave functions of the form 

s 



1 



1 + z 



(6) 



where f(z) is a polynomial of degree up to 2S + 1 in z. 

The differential area element on the surface of the sphere in terms of the stereographic 
coordinates x and y is 

AR 2 



dA 



:i + u 



2\2 



dxdy 



(7) 



and so the scalar product between any two polynomials / and g in this Hilbert space is 

4R 2 



if, 9) 



(1 + 1* 



2A2S+2 



f*g dx dy. 



(8) 



With this definition of the scalar product it is straightforward to derive the following identity 



by repeatedly integrating by parts, exploiting the fact that ^(f( z )*) — 0, 



/' dz n9 l ~ (2S + 2-n) 



9 ■ 



(9) 



i + kl 2 / / 

This result immediately implies the following spherical generalization of the z* — > d/dz rule 
of Girvin and Jach to the sphere, 
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U + (25 + 2)! cfe™ 1 ; 

Thus, to project any polynomial f(z, z* / (1 + \z\ 2 )) into the lowest Landau level one simply 
orders each term so that the z*/(l + l-^ 2 )^ all sit to the left, then replaces these factors by 
derivatives with respect to z according to (|10|). 

B. N particles and the fractional QHE 

The Hamiltonian for N spin polarized electrons on the sphere interacting via the Coulomb 
repulsion is then 

N 

H = J2Ti + V. (11) 

i=l 

The interaction energy V is 

v = -r£. + l2.-»* (12) 

k r 'i 2R R ' y > 

where is the chord distance between a given pair of electrons on the sphere. Here the 
sphere is assumed to have a uniform compensating charge density with total charge Q. When 
considering a homogeneous state the appropriate background charge is Q = Ne for which 

e 2 1 N 2 e 2 

V = Y—--—^. (13) 

h r « 2 R 

The spherical analog of the Laughlin state |§ at v — 1/q occurs when q(N — 1) = 25* 
and in the Wu-Yang gauge is 

s 



z,y. (M) 



The quasihole wave function corresponding to a single charge — e/q defect at the top of the 
sphere (z = 0) is 

k \ ^ I fc l / i<j 
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III. QUASIELECTRON WAVE FUNCTIONS 



A. Composite Fermion Quasielectron Wave Function 



According to Jain's composite fermion approach H a given fractional QHE wave function 
at filling fraction v = p/(kp + 1) where k = 2,4, •■■ is found by first constructing the 
corresponding integer QHE wave function at vcf — P and then multiplying by a Jastrow 
factor which ties k vortices to each electron. The state is then explicitly projected into 
the lowest Landau level. Denoting the Slater determinant corresponding to the effective 
integer QHE state, with the overall 11(1+ |z| 2 ) -5 factor removed, as <&cf, the corresponding 
fractional QHE states are 



(16) 



where Plll is the projection operator onto the lowest Landau level. 
For the Vcf — 1 ground state the Vandermonde determinant 



CF 



1 z\ z\ 



1 Zn zj 



„AT-1 



JV-1 



n( 

Kj 



(17) 



corresponding to one filled 'pseudo'-Landau level of composite fermions, gives for ip the 
Laughlin wave function ( p!4|) for v = 1/q where q = k + 1. If we remove a composite fermion 
from the lowest 'pseudo'-Landau level then 



CF 



Z\ z 1 • ■ ■ Z 1 



7,r 7 2 ■■■ 7 N ~ X 



n^n(^ - 



(18) 



k i<j 



and ip is identical to ( |T5| ) and describes a state in which a single quasihole sits at the top of 
the sphere. 

We now consider the quasielectron wave function constructed using this approach. If we 
introduce a composite fermion into the first excited 'pseudo'-Landau level then 



'>CF 



N-2 



1 Zi ... Z 1 



1 -y -y N - 2 
1 Z N ... Z N 



1 + 



l+kiv| 2 

and the corresponding physical electron wave function is 



(19) 



/,| IJ ^ HkJ 2 



n(^ - 
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1 *1 



1 Ztv 



JV-2 
*1 



i+kil 5 



JV-2 



(20) 



^ l+kivl 2 

This wave function can be projected into the lowest Landau level as follows. First rewrite 
tpqjf by pulling the Jastrow factor and the projection operator into the last column of the 
determinant, 



n r 



+ \ z k\ 



1 z l 



„AT-2 



„iV-2 



PlA,L 1+ \ Zl \2 Ili<j( z i z j) q 



1 Z N ... Z N ^LLLiq^p 



(21) 



Tli<j( z i z j) Q 

This can be done here because the cofactor associated with the n th element of the N th column 



does not contain z n . Thus, when we do the projection we need only project each element of 
the matrix separately. Following the procedure outlined in Sec. IIA this projection gives 



Thus ipqf can be rewritten, up to an irrelevant normalization constant, to give 



i<3 



2S + 2 dz x 



9 n^-^r^^ife-^E 1 
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(22) 
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z i L,v£i 
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(23) 



One can simplify things further by expanding the determinant in a cofactor expansion 
down the N th column. The cofactors are then all Vandermonde determinants and the final 
expression for the Jain quasielectron is 

1^1 



'qe 



En 



E 



n k^n Zk Zn l^n Zl Z ' a 



-1>; 



gs- 



(24) 



In this form 2 can be sampled by usual variational Monte Carlo techniques with each 
Monte Carlo step taking order iV instructions, rather than N 2 for a usual determinant. 



B. Laughlin Quasielectron Wave Function 

The generalization to the spherical geometry of the quasielectron wave function intro- 
duced by Laughlin is 

Straightforward Monte Carlo sampling of \ip qp \ 2 is not possible because of the explicit deriva- 
tives with respect to the electron coordinates. To compute the energy of this state we there- 
fore follow the procedure of Morf and Halperin |5||J which, for completeness, we review 
below. A more detailed discussion can be found in ||. 

Following we first take the absolute square of the wave function to obtain 

in/ ! V JLJL 1 1 r 

\6tn Af\ 



n(i7y>*)n>'- 2 ' 12 * (26) 

The expectation value of any operator O depending only on the coordinates {(xi,yi)} is 
then 

which, after integrating by parts twice in the numerator and the denominator can be rewrit- 
ten 

^ = JP(zi,...,z N )dYli dxjdyj 

fP{z 1 ,...,z N )U i dx i dyi ~ 



where 



i \ 2s+4 /. i 



and 

25+4 



The chord distance between any two points on the sphere is given by = R\z{ — 
Zj\/ \J(1 + |^| 2 )(1 + \zj\ 2 ) and so for the Coulomb interaction the operator O is 



e 2 Jl + lzAiJl + lzA 2 
Vcoui. = ^ _ V , ' J ' (31) 

It is then straightforward to compute O and evaluate (|28|) by usual variational Monte Carlo 
techniques. 



IV. RESULTS 

The excitation energies of isolated quasielectron and quasihole states have been obtained 
using the trial wave functions reviewed in Sees. II and III. Following || we have computed 
the proper energies, meaning that the relevant ground state energies are computed with 
monopole strength 2S = q(N — 1) and background charge Q = Ne, while the energy of the 
quasielectron (quasihole) excitations are obtained keeping R and iV fixed and decreasing 
(increasing) the monopole strength according to S — > S — 1/2 (S —>■ S + 1/2). In addition, 
following H| and 0, we have shifted the background charge when computing the quasielec- 
tron (quasihole) energies taking Q — > Q — e/q (Q — > Q + e/q), in order to compensate the 
charge density of the bulk of the wave function. This eliminates a finite size correction of 
±(e/q) 2 l/R ~ 0(~7=), a correction which was not included in and which may account 
for the slightly different results obtained here. Our results for the proper energies of the 
quasiholes, the Laughlin quasielectrons, and the composite fermion quasielectrons for filling 
fractions v = 1/3,1/5 and 1/7 are shown plotted vs. 1/N in Fig. 1. 

The quasihole energy gaps obtained here, A q ' h ', extrapolated to the iV — > oo limit are 
given in Table I, together with the v = 1/3 result of Morf and Halperin || and the ex- 
trapolated exact diagonalization results of Fano et al. ||. As stated above, the discrepancy 
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between our results and those of Morf and Halperin is most likely due to the 1 / y/~N finite 
size correction we have eliminated by modifying the background charge density. 

The quasielectron energy gaps obtained using both the Laughlin trial state, A q L e ' , and 
the composite fermion state, A q CF are givin in Table II. Again, in comparing the present 
v = 1/3 result for A q L e ' with those of || we note a slight discrepancy which we attribute 
to the 1/yN finite size correction we have included here. For u = 1/3,1/5, and 1/7 the 
composite fermion excitation energy A q CF is consistently about 10 % lower than the Laughlin 
excitation energy, as can be seen clearly in Fig. 1. Note that the composite fermion result is 
also in slightly better agreement with the extrapolated exact diagonalization results of Fano 
et al. |§. 

Finally, Table III gives results for the total energy gap for creating a well separated 
quasielectron quasihole pair, A = A 9 "' 1 ' + A 96 '. The results are again compared with those 
of U for the Laughlin energy gap as well as those of |7j for the composite fermion energy gap. 
For v = 1/3 our extrapolated energy gap using the Laughlin quasielectron is A F = 0.110 ± 
0.002, roughly 20 % larger than the earlier estimate of Morf and Halperin ||. Our improved 
calculation gives a Laughlin energy gap which is slightly larger than the corresponding 
energy gap computed using the composite fermion quasielectron, Aqf = 0.106 ± 0.002. 
This is consistent with the results of Girlich and Hellmund f| using the short-range V\ 
model; however, we find here that for the Coulomb interaction the energy gaps obtained 
using these two wave functions are essentially the same. Table III gives similar results for 
v = 1/5 and 1/7. We therefore conclude that that both the Laughlin and composite fermion 
quasielectron wave functions provide adequate descriptions of the physical quasielectron, 
though the energy gap obtained using the composite fermion quasielectron is slightly closer 
to extrapolated exact-diagonalization results for all the filling factors we have considered. 

To summarize, the quasielectron and quasihole excitation energies in the fractional QHE 
have been calculated for v = 1/3,1/5 and 1/7 by variational Monte Carlo. Results have 
been obtained using the quasielectron states originally proposed by Laughlin, as well as the 
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fully projected composite fermion quasielectron states proposed by Jain. We have improved 
on earlier estimates of the excitation energies of the Laughlin states at v = 1/3 in 
order to show that the composite-fermion energy gap is actually slightly lower than the 
Laughlin energy gap, consistent with the results of Girlich and Hellmund ||. Results for the 
energy gap using Laughlin's quasielectron for v = 1/5 and 1/7, obtained here for the first 
time, show that, as for v = 1/3, for Coulomb interactions the energy gaps obtained using 
the Laughlin and composite fermion quasielectron wave functions are essentially the same, 
though those obtained using the composite fermion quasielectrons are slightly smaller, and 
closer to extrapolated exact-diagonalization results of Fano et al. 0, than those obtained 
using the Laughlin quasielectron. 
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TABLES 



TABLE I. Quasihole energ 


y for the fractional QHE with v = 1/3, 1/5, and 1/7. The Monte 


Carlo results of Morf and Halperin || and extrapolated exact diagonalization results of Fano et al. 


H are given for comparison. 








V 


A q.h. 


A^\Ref.[\\) 


A* h iRef.[2}) 


1/3 


0.0279(12) 


0.0224(16) 


0.0264 


1/5 


0.0092(6) 




0.0071 


1/7 


0.0038(4) 






TABLE II. Quasielectron 


energy for the fractional QHE with v = 1/3, 1/5, 


and 1/7. The 


Monte Carlo results of Morf and Halperin Q and extrapolated exact diagonalization results of 


Fano et al. || are given for comparison. 










At c -(Ref. §) 


A^-(Ref. §) 


1/3 0.0779(10; 


1 0.0825(12) 


0.075(5) 


0.0772 


1/5 0.0166(6) 


0.0191(6) 




0.0173 


1/7 0.0063(4) 


0.0070(5) 






TABLE III. Total gap for the fractional QHE with v = 1/3, 1/5, and 1/7. The Monte Carlo 


results of Morf and Halperin |] 


J and Bonesteel |@|, 


and extrapolated exact diagonalization results 


of Fano et al. || are given for comparison. 






v A C f 


A L 


A L (Ref. §) A CF ( U) 


A(Ref. §) 


1/3 0.1058(16) 


0.1104(17) 


0.092(4) 0.106(3) 


0.1036(2) 


1/5 0.0258(9) 


0.0283(9) 


0.025(3) 


0.0244(3) 


1/7 0.0101(6) 


0.0108(6) 


0.011(3) 
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FIG. 1. Proper energies for the Laughlin quasielectron (solid circles), composite fermion quasi- 
electron (solid squares), and quasihole (solid diamonds) for filling fractions = 1/3, 1/5 and 1/7, 
plotted vs. 1/N. 
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